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Abstract 

Using density functional theory we calculate the density profiles of a binary solvent adsorbed 
around a pair of big solute particles. All species interact via repulsive Gaussian potentials. The 
solvent exhibits fluid-fluid phase separation and for thermodynamic states near to coexistence 
the big particles can be surrounded by a thick adsorbed 'wetting' film of the coexisting solvent 
phase. On reducing the separation between the two big particles we find there can be a 'bridging' 
transition as the wetting films join to form a fluid bridge. The potential between the two big 
particles becomes long ranged and strongly attractive in the bridged configuration. Within our 
mean-field treatment the bridging transition results in a discontinuity in the solvent mediated 
force. We demonstrate that accounting for the phenomenon of bridging requires the presence of a 
non-zero bridge function in the correlations between the solute particles when our model fluid is 
described within a full mixture theory based upon the Ornstein-Zernike equations. 

PACS numbers: 
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I. INTRODUCTION 



Big solute particles (e.g. colloids) immersed in a solvent of smaller particles interact 
with each other by an effective potential which is the sum of their direct interaction and a 
solvent mediated (SM) potential. Even when the direct interaction consists solely of two- 
body terms, the SM potential contains higher body contributions, of all orders which are 
determined formally by integrating out the solvent degrees of freedom. This conceptual 
framework yields, in principle, a much simpler effective Hamiltonian which involves only 
the coordinates of the big particles.— In certain systems the two-body term in the SM 
potential may dominate completely the corresponding direct interaction. A well-known 
example is a suspension of big hard-sphere colloids in a solvent of small hard-spheres. There 
the SM potential between the colloids is termed the depletion interaction and this is the 
only contribution to the effective potential for separations greater than the big hard-sphere 
diameter.— In the case of a (non-hard) solvent which is at a state point near to fluid-fluid 
phase separation, big solute particles can be surrounded by a thick adsorbed 'wetting' film 
of the coexisting solvent phase.— If two such big particles become sufficiently close, there 
can be a 'bridging transition' as the wetting films surrounding the two big particles join to 
form a fluid bridge of the wetting phase - see for example Ref. U and references therein. In 
wet granular media these bridging (or capillary) forces lead to strong and very short-ranged 
interactions. For small solutes these effective attractions become rather long-ranged with 
respect to the dimension of the solute as is known from tip-substrate interactions in atomic 
force microscopy— Long-ranged attractive interactions are also surmised for hydrophobic 
molecules in water at ambient conditions.— Bridging is also a possible mechanism for driving 
colloidal flocculation. 7 

In previous work,^ the wetting of a binary solvent around a single big particle and 
the influence of these thick adsorbed films on the effective SM potential between two big 
particles, was investigated for a particular model fluid, namely the generalisation to mixtures 
of the Gaussian core model (GCM) . 1 ! 8 ' 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ' 18 A Gaussian potential provides a 
good approximation for the effective potential between the centres of mass of polymers in 
solution . 1 i 19 i 20 The approach to calculating the SM potentials was based upon the theory 
developed by Roth et al£> - henceforth referred to as the 'insertion method'. The insertion 
method works within the framework of density functional theory (DFT)2I and uses as input 
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the density profiles calculated around a single big particle in order to calculate the SM 
potential between a pair of big particles.— Although the insertion method is formally exact, 
in practice one must employ an approximation for the free energy functional of the mixture 
of big and small particles. 2 For state points near to coexistence we found thick adsorbed films 
around the big particles resulting in long ranged, strongly attractive SM potentials whose 
range was determined by the thickness of the wetting film. However, using the insertion 
method, we were unable to detect any direct sign of bridging in the SM potential. 8 ' 9 

The present work can be viewed as going a significant step further than Refs. UJ Here 
we investigate the same system: two large solute Gaussian particles, immersed in a binary 
GCM solvent near to phase separation. However, whereas the previous work used the elegant 
insertion method, the present work can be viewed as the 'brute-force' approach to the 
problem. Using an accurate DFT for the binary GCM solvent of small particle o 1 ! 8 ! 9 ! 16 ' 17 we 
calculate explicitly the solvent density profiles around a fixed pair of the big GCM particles, 
treating the latter as external potentials, and determine the resulting grand potential. By 
repeating this calculation for a range of values of the separation between the centres of the 
two big particles we obtain the SM potential. We find, within the present (mean-field) DFT 
approach, that when thick adsorbed films are present there can be a bridging transition as 
the separation between the two big particles is decreased, i.e. the formation of a bridging 
configuration gives rise to a discontinuity in the derivative of the SM potential. We believe 
that this is the first time a non-local DFT has been used to calculate bridging density profiles 
and the resulting effective potential between two particles. Bridging has been investigated 
previously within coarse-grained, local DFT in the recent study of Stark et alM for big hard 
spherical colloids immersed in an isotropic liquid crystal host close to the isotropic-nematic 
phase boundary. Similarly, Andrienko et a/.— calculated bridging density profiles of a solvent 
adsorbed between a big colloid and a planar wall using a local DFT. 

We also investigate the SM potential between two big GCM particles in a region of 
the solvent phase diagram near the binodal but lying outside the single particle thin-thick 
adsorbed film transition line^ where a single big particle does not have a thick adsorbed 
'wetting' film of the coexisting solvent phase around it. Adsorption still influences strongly 
the SM potential. We find an analogue of capillary condensation; as the two big particles 
become sufficiently close, the composite object is large enough to induce condensation of the 
coexisting solvent phase around the pair of big particles. This effect is somewhat different 
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from that which can occur between two big hard-core particles in a solvent near to coexis- 
tence. When a pair of such particles are sufficiently close, a bridge of the coexisting phase 
can condense in the gap between the two big particles, without there being thick 'wetting' 
films adsorbed on each of the big particles.— In the present soft core system the strong ad- 
sorption is not confined to the space between the big particles, rather it extends through the 
whole region in which the two big particles are situated. This local condensation also results 
in a jump in the SM force between the two big particles with the SM potential becoming 
strongly attractive for small separations. 

In the final part of the present work we relate our results for the SM potential to an 
approach for calculating the SM potential based upon the mixture Ornstein-Zernike (OZ) 
equations . 26 i 27 i 28 By solving the OZ equations together with a closure relation one can cal- 
culate the various fluid correlation functions. It is well known that if one makes a diagram- 
matic expansion for the fluid correlation functions the hyper-netted chain (HNC) closure 
approximation neglects a certain class of (bridge) diagrams which, taken together, is termed 
the bridge function, 26 We show that in order to account for the phenomenon of bridging 
of solvent between big particles within an OZ approach to the fluid structure, one must 
incorporate an accurate theory for the bridge- diagrams. 

The paper is laid out as follows: In Sec. |Tl]we describe briefly our model fluid, the GCM, 
and the DFT used to calculate the solvent density profiles and the SM potential between 
two big solute GCM particles. Section ITTTI presents results for the density profiles and SM 
potentials in the regime where there are thick adsorbed films around a single big particle 
resulting in a bridging transition when two big particles are sufficiently close together. In 
Sec. HVI we present a simple analytic 'capillarity' approximation which describes qualitatively 
the onset of the bridging transitions that we find. Section IS describes the effect of the 
formation of a thick adsorbed film around a pair of big particles, in the portion of the phase 
diagram where there is no thick film around a single big particle and Sec. IVII describes our 
demonstration that bridging between big particles is related to the bridge-function. Finally, 
in Sec. IVII| we discuss our results and draw some conclusions. 
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II. MODEL FLUID AND SM POTENTIALS 



We determine the SM potential between two big (B) Gaussian particles immersed in a 
binary solvent of smaller Gaussian particles. The GCM, in which the particles of species i 
and j interact via purely repulsive Gaussian potentials 

%(r) = exp(-r 2 /-Ry)> (1) 

r-i 

is a simple model for polymers in solutio n 1 ' 12 ' 13 ' 16 ' 17 ' 20 (in particular, Ref. [l| provides a good 
general introduction to the GCM). For the binary GCM solvent we choose pair potential 
parameters corresponding to a binary mixture of polymers of length ratio 2:1, as were used in 
previous work on this model fluid . 8 ' 9 ' 16 ' 17 The values are R22/R11 — 0.665, R12/R11 = 0.849, 
/3en = /5e 2 2 = 2.0 (f3 = l/k B T) and e^/en = 0.944. R u is the basic length scale in 
the system. For this choice of parameters the binary mixture exhibits fluid-fluid phase 
separation. The phase diagram of this binary solvent is plotted in the total density, p° = 
Pi + P°>> ver sus concentration, x = pPj/p , plane in Fig. [T] (p® are the bulk densities of the 
small particles of species v — 1, 2) - see also Ref. [rJ 

The SM potential between two big particles, labelled A and B, with centres at and 
r B , separated by a distance h, is given by the difference in the grand potential: 

W AB (h) = Q(\r A - r B \ = h) - Q(\r A - r B \ = 00). (2) 

This result can be re-expressed (trivially) in terms of excess grand potentials, uj l ex = Q — fib, 
where i = A, B and fib is the grand potential of the bulk solvent in the situation where there 
are no big particles present. Then, 

W AB (h) = ui B (\r A - TB \ =h)- ui - ug. (3) 

Note that LU l ex , the excess grand potential for inserting a single big particle of species i, is 
equal to fi l ex , the excess chemical potential of big species i in the limit of the bulk density 
of this species p° — > O.^&S The effective pair potential between two identical big particles is 
then the sum of the bare interaction v BB (r) and the SM potential: 

v e B f i(h) = v BB (h) + W BB (h). (4) 

Recall also, that 

v BB( h ) = -k B T\ng BB (h), (5) 
5 
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FIG. 1: The bulk phase diagram for a binary mixture of GCM particles with eyzjtw = 0.944 
and R22/R11 = 0.665, equivalent to a mixture of two polymers with length ratio 2:1 (see also 
Ref.Q). p° is the total density and x is the concentration of the smaller species 2. The solid line 
whose ends are denoted by filled circles is the thin-thick adsorbed film transition of the binary fluid 
adsorbed around a single big GCM particle with pair potential parameters (3eBi = 1-0, (3es2 = 0.8, 
Rbi/Ru = 5.0 and Rbz/Rii = 4.972 - see Ref. ;9. It meets the binodal at the 'wetting point' 
(upper circle) with x = 0.975 and p°Rfi = 10.1 (note these values differ slightly from the result 
quoted in Ref . - see footnote^ ) and terminates at a critical point (lower circle) with x = 0.94 
and p Rfi = 7.5. The solid line whose ends are denoted by filled squares is the thin-thick adsorbed 
film transition of the binary fluid adsorbed around a composite pair of the same big GCM particles 
at zero separation h = 0. This transition line meets the binodal (upper square) at x = 0.995 and 
p°Rf 1 = 14 and terminates at a critical point (lower square) with x = 0.973 and p°Rf 1 = 8.5. 

where qbb is the big-big radial distribution function in the limit of the big particle bulk 
density p B — ► 0. In the present work we use DFT to obtain the quantities uj^{\ya — yb\ = h) 
and u l ex . 

In DFT one calculates the solvent one body density profiles, {p v (r)}, for a given set of 
external potentials, {V^(r)}, by minimising the grand potential functional:— 

where \x v are the chemical potentials for the two species, v = 1, 2, of solvent particles. The 
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ideal gas part of the intrinsic Helmholtz free energy functional is 

JwKM] = k B T^2 I dr p u (v) [ln(A 3 ^(r)) - 1], (7) 

where A„ is the thermal de Broglie wavelength of species u, and ^^[{p^}] is the excess 
part of the intrinsic Helmholtz free energy functional. Minimising (JHJ) together with (JJJ) one 
obtains the Euler-Lagrange equation 

= k B ThxKlp v {v) - k B Tc^(r) - fi v + K(r), (8) 

where 

c , 1)(r) = _^^ t (9) 

is the one body direct correlation function, which is a functional of {p v }- In an exact 
treatment the density profiles {p u } satisfying JHJ) would yield the exact grand potential Q 
as the minimum of fiy.— At this point we also recall that the two-body direct correlation 
functions are given by the second functional derivative^ 



5p v (r)5p e (r A 



For the GCM the following approximate excess Helmholtz free energy functional turns 
out, despite its simplicity, to be remarkably accurate at high densities p°Rf l > S : 1 ^ 12 ' 1 ?^ 2 ? 

F™ A [{pi}\ = \ E / dr / dr ' ^ W Pd*'M\r - (11) 
v£ J J 

where v„^(r) is the pair potential between the small solvent particles of species v and 
£, given by Eq. (0). The functional, Eq. (flip, is that which generates the RPA closure: 
c vl' RPA ( r : r ') = — ^fjf (l r — r 'D) f° r the P a i r direct correlation functions .^ 12 ^ 3 ' 16 The higher 
the density, the more accurate is the RPA for this soft core model.— 

In the present work we choose the external potential to correspond to two fixed big 
Gaussian particles of the same size, separated by a distance h: 

K(r) = e^exp(-(r + h/2) 2 / J RL) + e^exp(-(r-h/2) 2 /i?L), (12) 

with v = 1,2 and where h is a vector along the z-axis, with |h| = h, i.e. the centres of 
the big particles are at z = ±/i/2. Throughout the present study we choose the external 
potential parameters to be (3e B i = 1.0, (3e B 2 = 0.8, R B xJ Rn = 5.0 and Rm/Ru = 4.972, 



FIG. 2: The density profile of a one component fluid of Gaussian particles, with bulk density 
P\R\i = 6.9, around a pair of big Gaussian particles, whose centres are located on the z-axis a 
distance h/R\\ = 12 apart. The contours in the z — r plane correspond to pi(z, r) = 6.82, 6.84, 
6.86 and 6.88. 

the same values as those used for the big-small particle pair potentials in much of the 



work in Refs. With this external potential the solvent density profiles have cylindrical 
symmetry - i.e. the density profiles are functions p v {z, r), where the z-axis runs through the 
centres of the two big particles, and r is the radial distance from the z-axis. If the external 
potential on the solvent were exerted by hard particles, special care would be required to 
ensure that the hard boundary is compatible with the grid of the numerical calculations 
in order to avoid numerical artefacts in the contact density.— One would have to employ 
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either matching coordinate systems, such as the bispherical one used, e.g., in Refs. 
2J, or even more sophisticated finite-element methods with adaptive mesh-size.— One of the 



appealing features of the soft core GCM used in the present investigation is that we can avoid 
this problem and perform our calculations on a cartesian grid in cylindrical coordinates. 

In Fig. |21 we display a typical density profile for a one component solvent of particles of 
species 1, with the external potential given by Eq. (|12|) with h/R\i = 12. Having calculated 
the solvent density profiles for a given separation h of big particles, we can insert these into 
Eq. (0) to calculate the (excess) grand potential and the SM potential WebW from Eq. Q. 
In Fig. 01 we display the SM potential between two big GCM particles, calculated for a one 
component solvent with bulk density Pi-Rfi = 6.9, i.e. the state point corresponding to the 
profiles in Fig. |21 Figure El should be compared with Fig. 2 of Ref. 0. The open circles are 



the results from the present 'brute-force' calculation. The solid line is the result obtained 
using the insertion method, where one calculates only the solvent density profiles around an 
islolated, single big particle and then uses the general result^ 

PW BB (h) = 4\h - oo;p° B -> 0) - 4\h;p° B - 0), (13) 

i.e. one calculates the difference in the excess chemical potential between inserting the second 
big particle a distance h from the first and inserting it at h = oo. As emphasised in the 
Introduction, Eq. (fT3*|) is formally exact when we know the exact free energy functional for 
a mixture of big and small particles. Here we use the same RPA functional (JTTj) extended 
to include a third species B, in order to find an approximate c" B in Eq. (|T3*j) - see Ref. for 
more details. The results from the two different routes are almost indistinguishable for this 
point in the phase diagram, and generally for other state points where no thick adsorbed 
(wetting) films are present around the big particles. The dashed line in Fig. El corresponds 
to the analytic approximation for WsB(h) presented in Ref . 0: 

f3W p B u B re (h) = -(n/2f/ 2 Pe B1 p*R B1 ex P (-h 2 /2R 2 B1 ), (14) 

where p* = p\(3eBi/{\ + ^ 2 P^nRnPi) • The agreement between this approximation and 
the result of the full numerical DFT calculations is remarkably good. 



III. THE SM POTENTIAL WHEN THERE ARE THICK ADSORBED FILMS: 
BRIDGING 

We now consider the case when thick adsorbed films develop around the big GCM parti- 



cles. The circumstances in which this can occur are discussed in Refs. ISlJlJ In general there 
can be thick adsorbed films when the small solvent particles are in a state near to phase 
separation. For the present mixture, the big GCM particles favour species 1 of the small 
solvent particles, and so thick adsorbed films of the coexisting phase rich in species 1 can 
develop when the big particles are immersed in the solvent at a state point lying on the right 
hand side of the binodal, which is poor in species 1. In Refs. it was found that thick 
films develop via a thin-thick transition out of bulk coexistence. The locus of these transi- 
tions is shown as the solid line joining filled circles in Fig. ^ Note that this transition line 
meets the binodal at a 'wetting point' whose density is somewhat higher than that quoted 
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FIG. 3: The SM potential between two big GCM particles in a one component solvent of small 
GCM particles, with bulk density p\R\\ = 6.9. h is the separation between the two big particles. 
The solid line is the DFT insertion method results, the open circles are the results from the present 
'brute-force' calculation (the two are almost indistinguishable) and the dashed line is the analytic 
result, Eq. ((Til) , obtained in Ref. 0. 

in Ref. 0. This discrepancy is associated with the existence of metastable minima in the free 
energy.^ In Figs. 0] and |S] we display density profiles calculated for a pair of big particles 
immersed in a binary solvent of small GCM particles with bulk density p°Rfi = 8.5 and 
concentration x = 0.948, a state point near to coexistence, located inside the single particle 
thin-thick adsorbed film transition line (see Fig. Q). (Fig. 6 of Ref. M displays the solvent 
density profiles around a single big particle for this state point.) In Figs. 0] and El the centres 
of the big particles are a distance h/Ru = 17 apart and there are thick adsorbed wetting 
films around the big particles. However, in Fig. H] there is a fluid bridge between the two 
particles whereas in Fig. there is no fluid bridge. This second set of profiles corresponds 
to a metastable situation. For this state point the bridging transition occurs at a slightly 
larger separation h t /Rn = 17.4; this is where the bridged and unbridged configurations have 
equal grand potential. In Fig. El we display the SM potential WbbW for this state point. 
There are two distinct branches, corresponding to bridged and non-bridged configurations. 
For h > h t the unbridged configuration is the stable one, whereas for h < h t the bridged 
configuration becomes stable. Since the two branches of WBB{h) have different slopes there 
is a discontinuity in the SM force, — dWsBih) /dh, at h t , the separation where the transition 
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occurs. The extent of the metastable portions is substantial; these extend well beyond the 
equilibrium transition. This type of metastability, with accompanying hysteresis, was also 
found by Stark et alM in their recent study of the bridging of the nematic wetting film be- 
tween two colloids immersed in the isotropic phase of a liquid crystal. We display in Fig. [7| 
the SM potential calculated in the same way for a different point in the phase diagram, closer 
to the solvent bulk critical point, at a total density p°Rfi = 6.9 and concentration x = 0.88. 
This state point is also near to bulk coexistence (see Fig. [TJ. In both figures El and [7| we 
compare the SM potential calculated using the present 'brute-force' approach (solid lines) 
with the results obtained using the insertion method (dashed line) as described in Ref. |9j 
There is a significant difference between the results from the two methods; the insertion 
method does not capture the existence of two distinct branches of the grand potential. Thus 
it does not appear to include explicitly the effects of a bridging transition. The insertion 
method does predict very strongly attractive SM potentials, of a similar magnitude to those 
from full DFT, but does not yield the correct shape or range for Wrb(/i). In contrast we 
recall from Sec. |H] that in the regime where there are no thick adsorbed films, the results 
from the insertion method and the 'brute-force' method are in good agreement. 

IV. APPROXIMATE MODEL FOR THE SM POTENTIAL WHEN BRIDGING 
OCCURS 

In Ref. Owe found that when there was a thick adsorbed film around a single big particle, 
a good approximation for the excess grand potential of a single big GCM particle immersed 
in a binary GCM solvent of small particles is: 



where p c ° ex are the solvent bulk densities in the coexisting phase, i.e. the phase that forms 
the adsorbed film. I is the thickness of the adsorbed film (I ~ Rbv, but we determine its 
value by calculating explicitly via DFT, the density profiles around a single big particle) and 
7(Z) is the fluid-fluid surface tension, which we approximate by 7(00), the surface tension of 
the planar free interface (this is calculated using the approach presented in Ref. llfil ). The 
first term in Eq. (jl5j) is the excess grand potential for inserting a single big particle into the 
coexisting phase, obtained from the RPA bulk equation of state,— and the second term is the 
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FIG. 4: Density profiles p u (z,r), v = 1,2, for a solvent with total density p°Rfi = 8.5 and 
concentration x = 0.948, a state near to phase separation located inside the single particle thin- 
thick adsorbed film transition line (see Fig. The centres of the big particles are a distance 
h/Rn = 17 apart. Note the presence of thick adsorbed (wetting) films and the fluid bridge 
between the particles. The contours, plotted in the z — r plane, correspond to pi(z, r)/?^ = 1, 2 
and 3 and p2(z,r)Rf 1 = 2 to 8 in increments of 2. The bridged configuration is the stable one for 
this value of h/Ru. 

contribution from forming a spherical fluid-fluid interface. Generalising to two big particles 
we might therefore expect the following approximation to be reliable: 

2 

oog B (h) ~ 2J2^buR 3 Bu pT x + A(l,h) 7 '(l,h), (16) 

where A(l, h) is the surface area of the fluid-fluid interface between the adsorbed film of 
the phase rich in species 1 which develops around the two big particles and the bulk fluid 
rich in species 2. j'(l,h) is the surface tension, which we again approximate by 7(00), the 
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FIG. 5: Density profiles for the same state point and separation, h/R\\ = 17, as Fig. 21 but now 
there is no fluid bridge between the big particles. This configuration is metastable. 

planar fluid-fluid interfacial tension. A similar sharp-kink or capillarity approach was used 
in Ref. U to investigate bridging for very big hard-core solute particles that induce thick 
adsorbed (wetting) films but some new features arise for soft core systems. At first sight we 
might expect the first (single particle insertion) term in Eq. (Ilfij) to be inaccurate as h — > 0, 
when the big particles are strongly overlapping. However, this is not the case. When h = 
the first term in Eq. (|TB|) is accurate, since two big particles lying on top of each other result 
in an external potential that has the same form as that due to a single big particle with €b u 
twice the value for one of the big particles taken alone. In other words, if we take the first 
term in Eq. (|15j) and make the substitution e^ v — > 2es„, then we obtain the the same first 
term as in Eq. (jlfijl . Given this observation the first term in (fho|) should be accurate for 
both large h and for h = 0. Thus, by 'continuity' we expect it to be accurate for all values 
of h. The overall accuracy of Eq. (|TB|) should depend upon how accurately we determine 
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FIG. 6: The SM potential between two big GCM particles in a binary solvent of smaller particles for 
the same state point as in Figs. 0] and |SJ i.e. with total bulk density p°R\i = 8.5 and concentration 
x = 0.948. h is the separation between the centres of the two big particles. The dashed line is 
the result for Wbb{^) obtained using the insertion method, the dot-dashed line is the 'sharp- kink' 
result (see text, Sec. lIV|) and the solid lines denote the results from the 'brute-force' calculation. In 
the brute-force calculation, one finds that there are two branches for WsB{h) (see inset for more 
detail), each with a metastable portion. The branch with the smaller value of WsB{h) is stable. 
This corresponds to the configuration with no bridge for h > ht, and to the bridged configuration 
for h < ht- At ht/Rn = 17.4, where the two branches cross, there is a discontinuity in the gradient 
of Wbb(/i), i.e. there is a jump in the SM force at this separation. 

the surface area A(l, h) which appears in the second term. 

Using equations (fTBj). (fl3j) and we can obtain an expression for the SM potential: 

W BB (h) ~ [A(l,h) - 87r/ 2 ] 7 (oo). (17) 



We now present a simple model for A(l, h), (see also Ref. 32) which we expect to be reliable 
for values of h near to where the bridging transition occurs. 

When h ^> 21, i.e. no fluid bridge is present, then A(l, h) = 8tt1 2 , and Eq. (|T7|) gives 
WBB{h) = 0. When there is a bridge we approximate the end sections of A(l, h) by the 
surfaces of two sections of spheres with radius /, and the bridge surface by the surface 
generated by rotating the arc of a circle, of radius s, about the axis passing through the 
centres of the end sphere sections (the z-axis). We denote the width along the z-axis of the 
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FIG. 7: The SM potential between two big GCM particles in a binary solvent of smaller particles 
near to phase separation, with total bulk density p°Rfi = 6.9 and concentration x = 0.88. h is 
the separation between the two big particles. The dashed line is the result for the SM potential 
from the insertion method, the solid lines are the results from the 'brute-force' calculation and the 
dot-dashed line is the 'sharp-kink' result. In the inset we display a magnification of WsBih) for 
large h, showing the two branches crossing at ht/Ru = 22.7 and giving rise to a jump in the SM 
force. 

bridge section by 2w and the diameter of the bridge section at the mid point between the 
centres of the end sphere sections by 2d. The surface area of the two end spherical sections 
is 4ttZ(/ + h/2 — w) and the surface area of the bridge section is 4tts(s + d) arcsin(w/s) — 
Attsw. Requiring continuity of the surfaces where the end and bridge sections meet and also 
requiring continuity in the gradients at the point where these sections join, we eliminate s 
and d to obtain the following expression for the total surface area: 



We choose the value w = wq which minimises A, i.e. dA/dw\ w=Wo = and use this pre- 
scription for calculating A(l, h) with Eq. (JT7j) to calculate the SM potential between two 
big GCM particles at state points near to coexistence. For the case when the small particle 
solvent has a total density p°R\i = 8.5 and concentration x = 0.948, corresponding to the 
full DFT calculation of the SM potential in Fig. © we find that ^^7(00) = 0.830, and 




(18) 
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that the film thickness l/Ru ~ 7 (see Fig. 6 in Ref. Using these values in Eqs. (fTHJl 
and (fT7j) . we calculate the SM potential for this state point. The result is the dot-dashed 
line displayed in Fig. |U] which is in good qualitative agreement with our results from the 
full 'brute-force' calculation of the SM potential, particularly for values of h near to where 
the bridging transition occurs. We also used this simple approximation for the case when 
the small particle solvent has a total density p Q R\\ = 6.9 and concentration x = 0.88, cor- 
responding to the full DFT calculation of the SM potential in Fig. [7| For this state point 
f3R 2 u -f{oo) = 0.152, and l/R n = 9.6 (see Figs. 4 and 9 in Ref. M) and the SM potential 
is shown as the dot-dashed line in Fig. Again, the results are in qualitative agreement 
with those of the full calculation. In particular, this simple approach provides a surprisingly 
accurate means of estimating the value of h at which the bridging transition will occur. If 
we assume that bridging will only occur when WBB(h) < the resulting values of ht un- 
derestimate the results of the full calculation by only a few percent in both cases. Even for 
small values of h the results of the sharp-kink approximation for WBB{h) are of the correct 
magnitude. However, this approximation fails to reproduce the correct shape of WsBih) for 
small h. 



V. THICK ADSORBED FILMS ON COMPOSITE PARTICLES 

In the previous sections we considered only state points near the binodal where we know 
that a single big GCM particle is 'wet' by a thick adsorbed film of the coexisting phase 
rich in species 1, i.e. state points inside or below the single particle thin-thick adsorbed film 
transition line. However, there can also be pronounced effects on the SM potential due to 
the presence of thick adsorbed films for state points outside the single particle thin-thick 
adsorbed film transition line (see Fig.Q), where a single big particle immersed in the solvent 
does not develop a thick adsorbed film. When two big particles are sufficiently close together 
the resulting composite object can be sufficiently large that a thick film is adsorbed. This 
effect is somewhat analogous to the case for big hard-core solute particles, where for certain 
state points for which no thick adsorbed films are present, capillary condensation of the 
coexisting phase can occur in the space between the two big particles, provided these come 
sufficiently close together4^ We cannot strictly describe the phenomenon we observe as 
capillary condensation because the big particles that we consider in the present work have 
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soft cores. Nevertheless, the present phenomenon has a similar effect on the SM potential, 
i.e. there is a jump in the SM force on reducing the separation h. As mentioned above, 
this phenomenon occurs outside (but close to) the single particle thin-thick adsorbed film 
transition line. However, its occurrence is restricted to a particular region of the phase 
diagram. If one considers two big particles with full overlap (h = 0) one can calculate the 
thin-thick adsorbed film transition line for this composite object. This line is higher in total 
density than the corresponding single particle transition line (see Fig. ^) and serves as an 
upper bound for the regime where 'capillary condensation' occurs; the latter is restricted to 
the region betwen the two transition lines. 

The solvent density profiles around two big particles with h sufficiently small that this 
'condensation' has occurred are very similar in form to the profiles in Fig. 01 i.e. the 'conden- 
sation' does not just occur in the space between the two big particles, as would be the case 
with a pair of hard-core big particles. Rather, due to the soft core nature of the GCM, the 
adsorbed film spreads around the whole region in which the two big particles are situated. 

In Fig. IHlwe display the SM potential between two big GCM particles in a binary solvent 
of smaller particles with total bulk density p R\\ = 11 and concentration x = 0.983. This 
state point is located at bulk coexistence above the single particle thin-thick adsorbed film 
transition line but inside the transition line for the composite particle - see Fig. ^ For large 
values of h the SM potential calculated via the 'brute-force' approach is in good agreement 
with the results from the insertion method. At this state point the insertion method does 
not include any effects of thick adsorbed films since the inputs into this theory are the 
density profiles around a single big particle; for this state point a single big particle has 
no thick adsorbed film. However, as h is decreased the results of the full DFT calculation 
show that there is a discontinuity in the gradient of WsB{h) due to the formation of a thick 
adsorbed film around the two particles. The change in the SM potential is very pronounced; 
the potential becomes much more strongly attractive - see Fig. |H1 The insertion method 
(dashed line) accounts extremely well for the large h behaviour of the SM potential. It also 
describes accurately the metastable portion of Wbb()i) for h below the transition value. 
However, it fails completely to describe the stable, strongly attractive portion arising from 
the formation of the thick adsorbed film around the two particles; it underestimates the 
strength of the attraction by a factor of about 10. This is not too surprising given that this 
method inputs only the density profiles around a single big particle and that these exhibit 
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FIG. 8: The SM potential between two big GCM particles in a binary solvent of smaller particles 
with total bulk density p R\\ = 11 and concentration x = 0.983 (this state point is at bulk 
coexistence, outside the single particle thin-thick adsorbed film transition line, but inside the thin- 
thick adsorbed film transition line for a composite pair of completely overlapping big particles - see 
Fig-©- h is the separation between the centres of the two big GCM particles. The dashed line is the 
result for obtained using the insertion method and the solid lines denote the results from 

the 'brute-force' method. For this state point a single big particle does not develop thick adsorbed 
film, but when two big particles are sufficiently close together the resulting composite object can 
develop a thick adsorbed film. The two branches of Wbb{^) correspond to configurations without 
adsorbed films (stable at large h) and with films (stable at small h). These cross at h/R\\ = 10.2, 
resulting in a discontinuity in the gradient of WbbQt) and a jump in the SM force. 

no thick adsorbed films for this state point. 

VI. BRIDGING AND THE BRIDGE FUNCTION 

We recall that the SM potential Wbb(j') is related via Eqs. (jlj and (0) to the pair corre- 
lation function between solute particles, gBB{ r ), in a bulk ternary mixture which consists of 
a single, big solute species B and two solvent species, considered in the dilute limit of solute, 
Pb — > 0. Since integral equations are a standard tool to determine bulk pair correlation 
functions in the theory of classical liquids,-^ it is natural to analyse the SM potential within 
this framework. However, we recall from the outset that whilst integral equation theories 
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have achieved remarkable precision in the description of one-component bulk fluids, integral 
equation closure approximations are generally less reliable in multi-component mixtures, 
especially for situations where the size of one component becomes much larger than the 
others leading to the possibility of thick film adsorption or wetting phenomena or, in the 
case of hard-sphere mixtures, to depletion phenomena. 

Before we analyse Qbb{x) in the ternary mixture, it is instructive to point out some 
features of the two-component solvent mixture which provide a relationship between the 
hypernetted chain (HNC) integral equations and the RPA density functional used in the 
present work. Diagrammatic analysis yields the following standard relationships between 
the pair correlation functions in a homogeneous (bulk) mixture:— 

hij(r) - cg>(r) = J2 Pi f dr'Mk - r'\)c%(r% (19) 

k=l,2 ^ 

\n 9ij {r) + i3 Vij {r) = h tJ {r) - cg } (r) + ^(r), (20) 

where hij(r) = gij(r) — 1. The first equation is the OZ equation for binary mixtures, 
and the second provides the formally exact closure to the OZ equation in terms of the 
(generally unknown) bridge function bij(r). The bulk densities of the two solvent species 
are denoted by pi (k = 1,2). It is a special feature of the binary GCM (or related soft- 
core models) that its pair correlation functions are very well described within the HNC 
approximatio n 1 ' 12 ! 1 ^ 14 ' 1 ^ 1 ^ 1 ^ 1 ^ 2 ?!?? which amounts to setting bij(r) = 0. We denote the 
corresponding solution for the pair direct correlation function by c^' HNC (r). The relation to 
density functional theory follows by noting that the HNC equations, Eq. (j2Hj) with fey(r) = 0, 
are identical to the test particle equations obtained from a DFT with the excess free energy 
functional 

F™ C [{Pi\] = A ex ({p°}) + J2 ( drpr C A A (r) - 

i=l,2 J 

dr I dr' £ c^ HNC (\r - r'|)A Pl (r)Ap,(r'), (21) 

J ij=l,2 

corresponding to a Taylor expansion to quadratic order in Apj(r) about the bulk densi- 
ties. The test particle equations follow by choosing as external potentials the interparticle 
potential Vji(r), minimising the HNC grand potential functional with respect to Pj(r) and 
identifying g^r) = Pj(r)/pj. In Eq. (J2U), Api( r ) = pi(r) - p° and A ex ({p°}) denotes the 
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excess Helmholtz free energy of the bulk solvent. The HNC chemical potential is given by 

.HNC 



/3ft" 



E P ° f dr (Imomo - 4r w wi - c<?'" nc m) • (22) 

j=l,2 ^ V / 

Previous results for the GCM showed that the pair correlation functions obtained from the 
HNC were similar to those obtained from the RP AW^? and that the fluid-fluid binodals 
from the RPA and the HNC approximation were close to each other. 3 -- If one neglects the 
weak density dependence of cf^' HNC {r) and sets c^' HNC '(r) ~ c fi' RPA ( r ) — ~P v ji( r ) then 
one can show for the binary mixture, i = 1,2: 

- f* pa im] (23) 

where J r ^f A is the RPA functional defined in Eq. (JTTJ). 

We have seen in earlier sections that for a fixed big Gaussian particle exerting an external 
potential on the solvent close to coexistence, the RPA functional accounts for the formation 
of a thick adsorbed film. It also accounts for complete wetting at a planar wall. 17 Owing to 
the weak density dependence of c^' HNC (r) we also expect the HNC functional to describe 
thick film formation and complete wetting. (This is in sharp contrast to simple fluids of the 
Lennard- Jones type where the harshly repulsive core in the interatomic potential induces a 
strong density dependence of the direct correlation function c^(r) and the HNC functional 
(|2*T|) fails to account for complete wetting^* 3 ^). 

Explicit minimisation of the HNC functional for a binary mixture in the presence of 
an external potential due to a single solute particle yields the HNC solute-solvent integral 
equations. These can also be derived from the test particle equations of the HNC functional 
for the ternary mixture of binary solvent plus solute in the dilute limit of the solute, ps — > 0. 
This functional is linear in ps(r) and is, at most, quadratic in the other density profiles. It 
is given by 



^l C rn ~ + » H s NC j dr Mr) 

/ dr J dr' J2 ^ HN °^ ~ r'l)Api(r) p B (r') . (24) 



Here, PbHp^}) is the HNC insertion free energy (chemical potential) for inserting a single 
solute particle into the bulk solvent with densities p° (i = 1,2). Analogously to Eq. 
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H^NC j g gjy en \yy 



/dr (\h iB {r)[h iB {r) - c% HNC (r)) - c% HNC (r)) , (25) 

8=1,2 J \ / 



where the solute-solvent pair correlation function h iB (r) and the direct correlation function 

(2) HNC 

c iB ( r ) are determined by solving the solvent-solvent and solute-solvent HNC equations. 
In the dilute limit of solute, the solute-solvent direct correlation function satisfies the OZ 
equation 

cS(r) = h iB {r) I dr 'W|r - r'i)cf(r'), (26) 

.7=1,2 J 



for i = 1,2. In this treatment thick adsorbed films can develop around a big solute particle 
and this is manifest in the density profiles of the two solvent species and thus in h B j(r). It 
follows from Eq. (|26|) that information about thick films is fed into cj^(r). We can deduce 
that whenever thick film formation occurs, c^' HNC (r) can be very different from the RPA 
result —f3viB{r). 

We turn attention now to the solute-solute correlation functions. These are generated 
by employing ^Jtemi fi xrn g v BB( r ) as the external potential and minimising the grand 
potential functional with respect to ps(r). One finds 

\ng BB (r) + (3v BB {r) = £ p° I Av'h Bi {\v - v'\)cf^ HNC \r') . (27) 

i=l,2 ^ 

If one now employs the mixture OZ equations in the limit p B — > one obtains 

In g BB (r) + (3v BB {r) = h BB {r) - cf B {r) . (28) 

Note that the right hand side of Eq. (J27|) depends on the solute-solvent correlation functions 
h Bi (r) and and cf^ HNC \r). The former quantity is, essentially, the density profile of species 
i around a single big particle determined by minimising the HNC functional and the latter 
is given by the OZ equation ()26|) . One might expect both quantities to be given accurately 
by the HNC treatment. The g BB (r) resulting from Eq. (}2Tj) yields, via Eqs. and (jSJ), 
an SM potential which we refer to as W BB c {r) since this is consistent with the fact that 
9BB( r ) satisfies Eq. the HNC equation for big-big correlations; the latter sets the bridge 
function b BB (r) = 0. 

As the HNC inputs only the pair direct correlation functions c^- (r) of the small 
solvent species, which should be well described by their RPA counterparts, we adopt the 
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following procedure: determine the density profiles of the two small species around a single 
big particle by minimising the RPA grand potential functional, Eqs. © and ([lip, and 
use these as input for hsi(r), along with c^' RPA (r) for the solvent-solvent direct correlation 
functions, in Eq. (J2UJ) . The resulting cf B \r) are then used in Eq. (J27j) to calculate gBB( r ) and, 
hence, the SM potential - which should be very close to W^b ( r )- We find that W^ c (r), 
for large r, is almost identical to the branch of Wbb{^) obtained using the 'brute-force' 
DFT method presented in Sees. ITT1 and IHIl for which there is no fluid bridge. In other words, 
when there is no fluid bridge, i.e. for r = h > h t , where h t = h t (p\, p°) is the separation at 
which the bridging transition occurs, Wbb C ( r ) — Wbb{j") and we can infer that the HNC 
approximation bBB{ r ) — is valid. However, for h < h t we find W^ c {r) is very different 
from H^Bs(r), indicating that the bridge function obs^)-, omitted from this analysis, must 
be substantial for h < h t . Thus we have demonstrated that &bb(t) must play a significant 
role in determining the fluid structure when there is bridging. 

We conclude that the ternary HNC functional, Eq. ()24j) . describes correctly the bulk 
solvent-solvent correlations and captures thick film formation in the solute-solvent corre- 
lations with a vanishing solute-solvent bridge function, bis{r) = 0. For the solute-solute 
correlations when bridging is not present the HNC assumption &b_b(t) = remains accu- 
rate but this approximation fails completely when bridging is present. This means that a 
more sophisticated theory should include in the ternary functional terms proportional to 
Pb(i") A/?i(V) Apj(r") and higher orders. These will become important near the onset of the 
transition. 



VII. DISCUSSION AND CONCLUSIONS 



Using 'brute-force' DFT we have calculated the SM potential WsB{h) between a pair 
of big GCM particles in a binary solvent of smaller GCM particles. In particular, we have 
focused on the regime where the big particles are immersed in the binary solvent near to 
bulk phase separation, where thick films of the coexisting solvent phase adsorbed around the 
big particles influence strongly the SM potential. It is in this regime that we find bridging 
transitions. We show that the insertion method for calculating the SM potential used in 
Refs. 01^, and which is based on the ternary version of the RPA functional is unable 
to incorporate the effects of bridging. This method does provide an accurate approximation 
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for Wbb(/i) for solvent state points away from the binodal. The bridging that we find is 
of two types: i) that due to the joining of thick adsorbed films around the individual big 
particles, described in Sec. Illl[ and ii) that due to local condensation around a pair of 
particles, described in Sec. |V] Both result in a change in slope of Wbb{^) at a separation 
h = h t and therefore a jump in the SM force at h — h t . 

Within our mean-field theory, bridging manifests itself as a sharp (first-order) transition. 
However, this cannot be the case in reality since the bridging transition involves a finite 
number of particles and therefore fluctuation effects will round the transition (see discussion 

R 

in Ref. |4|). We can make a crude estimate of the extent of rounding effects by arguing that 
fluctuations should only be relevant when \W BB (h) — W BB hl '(h)\ < ksT, where W BB (h) 
denotes the branch of WsB{h) where there is a fluid bridge and W B ° B hT '(h) the branch without 
a fluid bridge. From this inequality we can obtain the width, 5h t , over which the transition 
at h t will be smeared. We find that Sh t /h t ~ 10 -2 for the state points corresponding to the 
SM potentials displayed in Figs. El and This measure of the rounding becomes smaller for 
solvent state points further removed from the bulk critical point. For bigger solute particles 
we also expect the extent of the rounding to become smaller. At first sight our estimate of 
the rounding may seem surprisingly small, bearing in mind that the size ratio between the 
big solute and small solvent particles is only about 7:1. However, due to the soft-core nature 
of the GCM fluid, the solvent density is high and the number of particles involved in the 
bridging transition can be large. This demonstrates one of the important differences between 
the soft-core GCM and more typical hard-core fluid systems: For hard-core particles one 
would not find thick adsorbed films of the solvent were the size ratio between the solute and 
solvent only 7:1. Typically the solute must be of order 50 or more times larger than the 
solvent particles for wetting phenomena to become significant - see also the discussion in 
Ref.El 

Our analysis in Sec. IVII demonstrates that in order to incorporate bridging into a full 
(ternary) mixture theory, one must implement an accurate theory for the fluid bridge func- 
tions; in particular for the solute-solvent and solute-solute bridge functions &B;(r) and 6bb(^)- 
b BB (r), the solute-solute bridge function, remains little understood but must play a crucial 
role when there are thick adsorbed films surrounding the big particles. That the bridge 
functions are required highlights the essential many-body nature of the effective interac- 
tion between the big solute particles. Hence, it is not surprising that the insertion method 
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combined with the ternary version of the RPA functional is unable to incorporate the 
effects of bridging on the SM potential. We re-iterate that the insertion method is formally 
exact; it is its use with an approximate functional which leads to neglect of the key features 
of bridging. In order to obtain insight as to what is required in a theory for the full mixture 
Helmholtz free energy functional which incorporates the effect of bridging, we consider the 
exact inhomogeneous Kirkwood-Hill formula*^ (recall Eq. Q): 



for the one-body direct correlation function of the big solute particles in the limit ps — > 0. 
VBv{r) are the big-small pair potentials and the parameter A, with < A < 1, is used to 'turn 
on' the effect of the inserted big particle via the potential \vbv{t). One calculates the solvent 
response through the inhomogeneous big-small pair distribution function gBv(?, r'; A), as A is 
increased from to 1. Combining Eq. (|29p with Eq. (fTHj) one obtains an exact expression for 
W B b{t), given by Eq. (70) of Ref. I| Consider the case when the solvent is near coexistence 
at a state point below the single big particle thin-thick adsorbed film transition line. If one 
calculates WbbW via Eq. (J2H1), then g Bu (r,r';\ = 0) will correspond to the distribution 
arising from a fixed single big particle located at r = — h/2 exerting an external potential on 
the solvent. This big particle will be surrounded by a thick adsorbed film. Then, 'turning 
on' the effect of the second big particle (by increasing A from zero) located at r = +h/2 
one could perhaps envisage the situation where there might be two 'jumps' in (^(r, r'; A) 
for a particular value of h = |h|. The first would be at A = Ai, when the potential \\Vbu{t) 
becomes sufficiently strong to induce condensation of the coexisting solvent phase around 
this second big particle. This 'jump' in g B v{^, r'; A) could then be followed by a second jump 
at A = A2 (Ai < A2 < 1), when a fluid bridge forms between the two big particles. That such 
complex phenomena must be described by Cg (r), which is obtained by taking one functional 
derivative of the excess Helmholtz free energy functional Eq. (|9*|). attests to the degree of 
sophistication required in the ternary mixture functional ^^[{pi}] if this is to incorporate 
bridging. By employing the RPA functional in the insertion approach one incorporates the 
effect of the thick film surrounding the first big particle via p„(r) in Eq. (j2HJ), but neglects 
the effect of the thick film around the second big particle by setting gBu{r, r'; A) = 1 for all 
A.— One is thereby unable to incorporate the effect of bridging on the SM potential. This 
point is highlighted further by the case described in Sec. El where a single big particle has no 
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thick adsorbed film and the bridging arises from condensation around a pair of big particles. 
In this situation all the information about bridging/wetting must be generated in c^\r) 
from a source other than the solvent density profiles around a single big particle, i.e. from 
subtle correlations in the inhomogeneous solvent. Incorporating such correlations is a tall 
order for a theory! 

The simple capillarity (or sharp-kink) approximation used in Sec. HVI to provide an ap- 
proximate theory for when bridging occurs between two thick adsorbed films seems to be 
quite good. The simple form, Eq. (j!7|) . taken with Eq. (fTHj) . is surprisingly reliable in deter- 
mining an approximation for ht, the separation between the big particles at which bridging 
occurs, as well as providing a reasonable approximation for the slope of WBB{h) near the 
onset of the bridged configuration, i.e. it provides quite a good approximation for the SM 
force at h ~ h t . The capillarity approximation is not reliable for small h. Here the shape of 
the SM potential determined from the capillarity approximation is completely wrong, and 
therefore the SM force obtained from this approximation will be completely unreliable - see 
Figs. IH1 and The 'brute-force' calculation shows that as h —>■ 0, the SM force — > 0, whereas 
the capillarity approximation shows the SM force tending to a non-zero constant value as 
h — > 0. We believe the origin of this failure lies in our simple approximation (jl8J) for the 
area of the fluid-fluid interface. 

One issue we have not raised so far is what does one take for vbb{ t ), the bare big-big 
pair potential? This does not enter our calculation of the SM potential, since it is only 
VBu(r), the big-small pair potentials, that are involved; the big particles are treated as 
external potentials. Therefore, in principle, f_Bs(r) could take any form, although choosing 
a bare potential with a hard core would be inconsistent with the soft-core nature of vbu{t). A 
Gaussian potential of the form given in Eq. would seem a natural choice for vbb{t). When 
one considers the GCM to be a simple model for polymers in solution, then the following 
empirical rules for the pair potential parameters apply between unlike species: ^ 1 ^ 1 ^ 2 ^ 1 Rf^j = 
+ R]j)/2 and e^j < en ~ e^. Therefore, the choices Rbb/Ru = 7 and €bb = ZksT 
would be consistent with the parameters we have used for the big-small pair potentials.^ If 
we employ a bare Gaussian potential with these parameters the big-big repulsion is negligible 
when compared to the attractive Wbb^)-, particularly when there are thick adsorbed wetting 
films present around the big particles. Thus, the resulting effective pair potential, v e J^(r), 
given by Eq. (@J) can be very strongly attractive. 



25 



If one were seeking to investigate experimentally the effects of thick adsorbed films and 
bridging between colloidal particles, one approach is to perform light scatteringexperiments 
in order to measure the second virial coefficient, B 2 .— As pointed out in Ref. y, B 2 , which 
measures the integral of — r 2 (exp[— /3v e J^(r)] — 1) should be very large and negative when 
adsorbed films are present. A rapid change to large negative values of B 2 upon changing 
the solvent state point should indicate the development of thick adsorbed films around the 
colloids, thereby influencing the SM potential. 9 Whether B 2 does show a rapid variation 
with composition in the neighbourhood of the thin-thick transition lines remains to be 
ascertained. 

Finally we note that, just as a pair of big particles with h = exhibits a thin-thick 
adsorbed film transition line, at total densities higher than the single particle thin-thick 
adsorbed transition line (see Fig.[TJ, there should also be a thin-thick adsorbed film transition 
line at even higher densities for three big particles whose centres coincide. This will have 
implications for the three-body interactions between the big particles. Furthermore, there 
may be other transition lines corresponding to four, five or more big particles completely 
overlapping. 
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